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Numerical range of a Hermitian operator X is defined as the set of all possible 
expectation values of this observable among a normalized quantum state. We ana- 
lyze a modification of this definition in which the expectation value is taken among 
a certain subset of the set of all quantum states. One considers for instance the 
set of real states, the set of product states, separable states, or the set of maxi- 
mally entangled states. We show exemplary applications of these algebraic tools in 
the theory of quantum information: analysis of fc-positive maps and entanglement 
witnesses, as well as study of the minimal output entropy of a quantum channel. 
Product numerical range of a unitary operator is used to solve the problem of local 
distinguishability of a family of two unitary gates. 

PACS numbers: 03.65.-w, 03.67.-a, 02.10.Yn 
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I. INTRODUCTION 

Expectation value of a Hermitian observable X among a given pure state belongs 
to the basic notions of quantum theory. It is easy to see that the set A of all possible 
expectation values of a given operator X among all normalized states forms a close interval 
between the smallest and the largest eigenvalue, A(X) = [\ m m, A max ]- 

In the theory of matrices and operators one calls such a set numerical range or field of 
values of an operator X, which in general needs not to be Hermitian [TJ |2]. Properties 
of numerical range are intensively studied in the mathematical literature [31 0], several 
generalizations of this notion were investigated [SHE], and its usefulness in quantum theory 
has been emphasized [9]. 

Let us introduce the set Q of all density matrices of size N, which are Hermitian, positive 
and normalized, Q := {p : p 1 = p > 0, Trp = 1}. If a given state is pure, p = |^)(^|, 
the expectation value reads TrpX = (i/j\X\i/j). Any density matrix can be represented as a 
convex combination of pure states. Hence for any operator the sets of its expectation values 
among pure states and among mixed states are equal. 

More formally, let X be an arbitrary operator acting on an A^-dimensional Hilbert space 
"Hat. Its numerical range, can be defined as 

A(X) = {TrXp, p G Q}. (1) 

The related concept of numerical radius 

r(X) = {\z\, zeA(X)} (2) 

is also a frequent subject of study 0, H] (c/. Table |ll| in Section |V|) . 

In this paper we analyze a modification of the standard definitions and ([2]). For any 
operator X, one defines its restricted numerical range 

A R (X) = {TrXp, P eQ R c Q}. (3) 

and the restricted numerical radius 

r R (X) = {\z\, zeA r (X)}. (4) 

The symbol fl R denotes an arbitrary subset of the set Q of all normalized density matrices 
of size N. Thus the above definition of the restricted numerical range is more general than 
the one studied in [21 [10] , in which a subset of the set of pure states was used. 
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Some examples of restricted numerical ranges are listed in Tab. |T} The range restricted 
to real states was recently discussed by Holbrook [TTJ, while the Liouville numerical range, 
in which the pure states of size M 2 reshaped into a square matrix form a legitimate density 
operator was analyzed by Silva [12]. The numerical range of a density matrix p restricted 
to the SU(2) coherent states gives the set of values taken by its Husimi representation 
- see e.g. [13]. Examples of restricted numerical radii can be found in Tab. [n] at the 
end of the paper. An very important example is the product numerical radius r®(X) = 
max{\{x(ji <g) . . . (g) tp m \X\ipi ® . . . <g) ip m )\ : G "H ni }, which coincides for m = 2 and X 
normal with the Schmidt operator norm ||-X"||gm introduced by Johnston and Kribs [H] . 

If the dimension of the Hilbert space is a composite number, N = KM, the space can be 
endowed with a tensor product structure 

Wn = Hk ® %m- (5) 

From a physical perspective this corresponds to distinguishing two subsystems in the entire 
system. One defines then the set of separable pure states, i.e. the states with the product 
structure, = \i/ja) ® \4>b)- 

Substituting this set into definition ^ of the restriced product range one arrives at the 
notion of product numerical range of an operator X, 

A®(X) = {(ip A <g> ip B \x\<ip A ® <M : \iPa) e H K , Mb) e Um) , (6) 

where both states \ipA) G Hk an d \iPb) G %m are normalized. 

The product numerical range can also be considered as a particular case of the decompos- 
able numerical range E] defined for operators acting on a tensor product Hilbert space. 
This notion was recently analyzed in [ToTfrS] . where the name local numerical range was 
used. In physics context the word 'local' refers to local action, so the unitary matrix with 
a tensor product structure, U(M) <8> U(K), is said to act 'locally' on both subsystems. To 
be consistent with the mathematical terminology we will use here the name "product nu- 
merical range", although a longer version "local product numerical range" would be even 
more accurate. Note that one may also use other restricted sets of quantum states as these 
mentioned in Table [H 

The main aim of this work is to demonstrate usefulness of the restricted numerical range 
for various problems of the theory of quantum information. This paper is organized as 
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TABLE I. Examples of restricted numerical range (NR): Ar(X) = {TiXp, p 6 ^_r}, where $1^ C SI 
denotes a subset of the set of all quantum states of size N. All pure states are assumed to be 
normalized, = 1, while all coefficients in the sums are non-negative. In each case eq. Q 

provides example of the corresponding notion of restricted numerical radius vr. 



follows. In Section [II] we review some basic features of product numerical range and present 
some examples obtained for Hermitian and non-Hermitian operators. Although we mostly 
discuss the simplest case of a two-fold tensor product structure, which describes the physical 
case of a bi-partite system, we analyze also operators representing the multi-partite systems. 



In Section III we study the notion of separable numerical range and other restricted numerical 
ranges of an operator acting on a composed Hilbert space. 



Key results of this work are presented in Sec. [TV] in which some applications in the theory 
of quantum information are presented. In particular, by analyzing a family of one-qubit maps 
we find the conditions under which the map is positive and establish a link between product 
numerical range of a Hermitian operator and the minimum output entropy of a quantum 
channel. The problem of fc-positivity of a quantum map is shown to be connected with 
properties of the numerical range of the corresponding Choi matrix restricted to the set VL^ 1 
of states with the Schmidt number not larger than k |19j . For k = 2, we point out that 
the question of distillability of an entangled quantum state is related to the numerical range 
restricted to the set fl^ 2 '. 

Furthermore, properties of product numerical range of non-Hermitian operators are used 
to solve the problem of local distinguishability for a family of two-qubit gates. In section 
|V| we present some concluding remarks and discuss further possibilities of generalizations 
of numerical range which could be useful in quantum theory. Proofs of certain lemmas are 
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relegated to the Appendix. 

II. PRODUCT NUMERICAL RANGE 

Quantum information theory deals with composite quantum systems which can be de- 
scribed in a complex Hilbert space with a tensor product structure [20]. When analyzing 
properties of operators acting on the composed Hilbert space ([5]), it is physically justified to 
distinguish product properties, which reflect the structure of the Hilbert space. 

If the physical system is isolated from the environment, its dynamics in time can be 
described by a unitary evolution l^') = U\ip) where U is unitary, UW — Ijv- In the case of 
a bipartite system, N = KM, one distinguishes a class of local dynamics, which take place 
independently in both physical subsystems, so that U = Ua®Ub, where XJ a G U(K), while 
Ub G U(M). From a group-theoretical perspective, one distinguishes the direct product 
U(K) x U(M), which forms a proper subgroup of U(KM). 

It is important to know which tasks, such as the discrimination of pure quantum states, 
can be completed with the use of local operations and classical communication. For this 
purpose, it is convenient to work with the notion of the product numerical range of an 
operator defined by Eq. (J6|. This algebraic tool can be considered as a natural generalization 
of the standard numerical range for operators acting on a tensor product Hilbert space. 

Note that the definition of product numerical range is not unitarily invariant, but implic- 
itly depends on the particular decomposition of the Hilbert space. This notion may also be 
considered as a numerical range relative to the proper subgroup U (K) x U (M) of the full 
unitary group U(KM). It is worth mentioning that product numerical range differs from 
so-called quadratic numerical range, also defined for operators acting on a composite Hilbert 
space [2T] . 

Consider the following problems arising in the theory of quantum information. 

• Verify if a given map $ acting on the set of quantum states is positive: Is > 
for all p > 0? 

• For a given observable X, defined for a bipartite system, find the largest (the smallest) 
expectation value among pure product states: What is max|0 Ai9iB ) (<pA, (Pb\^\<Pa, 4>b)7 
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• Check if two unitary gates U\ and U 2 acting on a bipartite systems are distinguishable. 
This is the case if there exists a product state \ip) = \4>Ai 0b) such that the states Ui\ip) 
and U 2 \ip) are orthogonal. 

• For a pair of two bipartite states a and p maximize their fidelity or the trace Trpcx by 
the means of local operations. 

This list of questions, of different difficulty levels, could be easily extended. All these 
problems have one thing in common: they could be directly solved, if we had an efficient 
algorithm to compute the product numerical range of an operator. Although in this work 
we are not in a position to go so far, we aim to show usefulness of this notion and present 
some partial results. 

A. Basic properties 

In this section we review some basic properties of product numerical range. Some of them 
were discussed by Dirr et al. [TS], while some other were established in [22J. 

For any operator X acting on a Hilbert space Hn, its product numerical range ^ forms 
a nonempty, connected set in the complex plane. However, this set needs not to be convex 
nor simply connected. Further properties of product numerical range include 

a) Subadditivity, A® ( A + B) c A® ( A) + A® (B) , 

b) Translation: for any a G C one has A (A + al) = A® (A) + a, 

c) Scalar multiplication: for any aGC one has A® (a A) = aA®(A), 

d) Product unitary invariance: A® (([/ ® V)A(U <g) V) ] ) = A® (A), 

e) If A is normal, then numerical range of its tensor product with an arbitrary operator B 
coincides with the convex hull of the product numerical range, A(A ® B) — Co(A®(A <g> B)) 

f ) Product numerical range of any A contains the barycenter of the spectrum, -A^ trA G 
A® (A). 

To analyze product numerical range of the Kronecker product it is convenient to make 
use of the geometric algebra of complex sets [23] . For any two sets Z x and Z 2 on the complex 
plane, one defines their Minkowski product, 

Z X M Z 2 = {z : z = z x z 2 , zi G Z x , z 2 G Z 2 } . (7) 



7 



Observe that this operation is not denoted by the standard symbol <g> in order to avoid the 
risk of confusion with the tensor product of operators. The above definition allows us to 
express the product numerical range of the Kronecker product of arbitrary two operators as 
a Minkowski product of the numerical ranges of both factors [2] , 

A (A® B) = A(A) H A(B), (8) 

This property can be directly generalized to an arbitrary number of factors. Thus the 
problem of finding the product numerical range of a tensor product can be reduced to 
finding the Minkowski product [231 121] of two or more numerical ranges. 

B. Hermitian case 

In the case of a Hermitian operator X = X* acting on Hn its spectrum belongs to the 
real axis. Labeling the eigenvalues in a weakly increasing order, Ai < A2 < ■ ■ ■ < Ajy, one 
can write the numerical range as an interval, A(X) = [Ai, Ajv], see e.g. [TJ. 

Let us assume, the Hilbert space has a product structure, Hn — Hk <8 'Km-, which 
implies a notion of a pure product state. Define the points A® in and A® ax as the maximal 
and the minimal expectation values of X among all product pure states. Then the product 
numerical range is given by a closed interval, A®pT) = [A® in , A® ax ]. If the spectrum of X is 
not degenerated to a single point, (which is the case iff X is proportional to identity), then 
A® in 7^ A® ax , so the product numerical range has a non-zero volume [22J. 

Making use of the lemma about the dimensionality of subspaces belonging to a composed 
Hilbert space of size N = KM which contain at least one separable state [25], one can get 
the following bounds for the edges of the product numerical range 

Amin — \k-x)(m-x)+i and A® ax > Xk+m-x- (9) 

These bounds, proved in [22], imply that in the simplest case of a 2 x 2 system (N = 4), the 
product numerical range contains the central segment of the spectrum, 

A(X) = [A 1; A 4 ] d A*{X) = K in , A^ ax ] D [A 2 , A 3 ]. (10) 

Similarly, for any Hermitian X acting on a 2 x K space the central segment of the spectrum 
[Xk, Xk+x] belongs to A®(X). In the case of a 3 x 3 system (N = 9), the product numerical 
range of X contains its central eigenvalue, A5 G A®(X). 
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1. Exemplary Hermitian matrix of order four 



Not being able to construct an algorithm to obtain product numerical range for an ar- 
bitrary Hermitian operator we shall study some concrete examples. Consider first positive 
numbers t, s > and a family of Hermitian matrices of order four 

f 2 t \ 
1s 
s -1 
\t -2 ) 



X, 



t,s 



(11) 



with the spectrum 



Then we can write 



{-v^TT, v^TT, V^Tl] . 



(12) 



(x\ <g> (y\X t)8 \x) <g> \y) = 2|^i | 2 |2/! | 2 + N 2 |y 2 | 2 - |x 2 | 2 |yi| 2 - 2\x 2 \ 2 \y 2 \ 2 

+2tRe[x* 2 x iy * 2yi ] + 2sRe[x* 1 x 2 y* 1 y 2 }. (13) 

Because in the case of Hermitian matrices the product numerical range forms a closed 
interval, we only need to find the upper and the lower bounds for the above expression. We 
have 

(x\ ® (y\X tjS \x) ® \y) < 2\ Xl \ 2 \ Vl \ 2 + |xi| 2 |y 2 | 2 - \x 2 \ 2 \yi\ 2 - 2\x 2 \ 2 \y 2 \ 2 

+2(t + s)|x 2 ||a;i||s/ 2 ||j/i|. (14) 

Because \xi\ 2 + \x 2 \ 2 = 1, we can put p = \x\\ 2 and q = \y\\ 2 with p, q ^ 0. This gives us 

(x\ ® ® \y) < 2pq + p{l-q) - (1 - p)q - 2(1 - p)(l - q) 

+2(t + s) s /p{l-p)q{l-q). (15) 

We want to maximize the above expression under the following constraints: < p < 1 and 
< q < 1. 

First we analyze the edge. On the edge (one of the variables p, q is or 1) the square 
root vanishes, the remaining part is convex and thus the extreme points are (p, q) G 
{(0,0), (0, 1), (1,0), (1, 1)}. Thus the maximum value on the edge is 2. If we assume that 
p, q {0,1}, we have to find zeros of appropriate derivatives. The extremum value is 
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A®(X , 5 ) J>^ 








A min 



1 



2 3 4 5 

s 



FIG. 1. Numerical range and product numerical range for matrices Xq s , which belong to the 



family ( 11 ). 



\/t 4 + lOt 2 + 9/2t for t > \/3. The lower estimate is obtained similarly. Thus the exact 
formula for the product numerical range reads: 



A®(X M ) = [-/(t + a),/(t + a)], 



(16) 



where 



/(*) 



(17) 



2 for t G [0, v^3) 

Vt 4 + 10t 2 + 9/2t for t e [y/3, oo) 
Note that the product numerical range depends only on the sum of the parameters s and 
t, whereas the numerical range depends on the values of both of them. The minimum and 
the maximum values in the numerical range and the product numerical range of the matrix 
Xt >a are compared in Fig. [TJ 

Let us consider a more general family of matrices for t, s > 

(a t \ 
b s 
s c 
\ t d J 

For given a,b,c, d, one can obtain a similar result as above, but in general the formulas are 
very complex due to the higher number of parameters. However, it is easy to obtain the 
following bound 

\f(8 + t),g(s + t)] CA%), (19) 



Y, 



t.s 



(18) 
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where 



and 



f(t) = min{min(a, b, c, d), -tvY , - ^*}> 



g(t) = max{max(a, b, c, d), -trlo.o + t^}- 



(20) 



(21) 



2. A tridiagonal Hermitian matrix 



Consider another family of Hermitian matrices of size four, written in the standard prod- 
uct basis, 



D 



2 

a* 



a 
\ b 



b* 1 c 



(22) 



where a and b are arbitrary complex numbers and c = xa for some arbitrary real number x. 

This family was introduced in j2S] as a useful example for studying block positivity. Here 
we deal with the product numerical range of D, but the two concepts are closely related, 
since a Hermitian matrix acting on a bipartite Hilbert space is block positive iff its product 
numerical range belongs to M + . Following the lines of [26], with some additional effort, one 



obtains an explicit result 



where 



A®(D) 



G 



' 2 



G 



G = \ (\a + c\ + v/|a-c| 2 + |6| 2 ) . 



(23) 



(24) 



3. Family of isospectral Hermitian operators 

It is instructive to study product numerical range for a family of Hermitian operators with 
a fixed spectrum and varying eigenvectors. Any unitary 4x4 matrix U may be represented 
in a canonical form: 

U = (V A ® V B ) U d (W A g» W B ) , (25) 
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where Va,Vb,Wa,Wb € U(2), while Ud is a unitary matrix of size four expressed in the 
form m\ 



U d {oti, 012, 013) = exp(i ^2 a k a k <g> a k ). 



(26) 



k=l 



Here a k denotes the Pauli matrices, and the three real parameters aj belong to the interval 

[Ml- 




0.4 



FIG. 2. Product eigenvalues and product numerical range (gray region) of the one-parameter (ai) 



family of matrices given by Eq. (27) with eigenvalues Ai = 0, A2 = 1/6, A3 = 1/3, A4 = 1/2. 



Consider a density matrix obtained from the diagonal matrix E(x\,X2, 23) = diag(xi, 22, £3, 1- 
X\ — X2 — X3) by a non-local unitary rotation, 



p(a x ,Xi) = U d EUl, 



(27) 



with c*2 = ct3 = 0. Figure [2] presents the dependence of its product numerical range as a 
function of the non-locality phase a±. 



4- Random Hermitian matrices of order four 

As shown in the above examples, the lower edge of the product numerical range of a 
Hermitian matrix X of order four is interlaced between its two smallest eigenvalues, A® in G 
[Ai, A2]. We have already seen that these bounds can be saturated, so the exact position of 
A® in is X dependent. However, following the statistical approach, one may pose the question 
how the edge is located with respect to both eigenvalues for a random Hermitian operator. 
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To analyze this problem we generated numerically a 5 x 10 5 random Hermitian matrices 
according to the flat (Hilbert-Schmidt) measure in the set Q of normalized density matrices 
of size N = 4. The joint probability distribution for the eigenvalues reads 

1*1 4 4 

P(A 1 ,A 2) A 3 ,A 4 ) = ^^l-^A^n^-A,) 2 . (28) 

j=l i<j 



3456 



By construction, the eigenvalues sum to unity, and this normalization sets the scale. It 
is possible to integrate out of the above formula any chosen three eigenvalues and obtain 
an explicit probability distribution for the last one. For instance the distribution for the 
smallest eigenvalue has the form 



P(Ai) = 60(1 - 4A X ) 14 6(Ai) 6(1/4 -Ai) 



(29) 



where Q(x) is the Heaviside function. 



P(X) 



A, 



.... 

A 2 



X 



i.o A 



FIG. 3. Probability density of eigenvalues Ax,A2,A3,A4 (dashed/dotted lines) and product values 



A 



min' max 



(dark histograms) for a random two-qubit density matrix, generated according to the 



Hilbert-Schmidt measure (28) 



Figure [3] presents the probability distributions for ordered eigenvalues, Ai < A2 < A3 < A4, 



obtained analytically by integration of (28). These distributions are compared with the 



distributions -P(A® in ) and -P(A® ax ) obtained numerically. As follows from (10) A® in is located 
between the two smallest eigenvalues, while A® ax is interlaced by the two largest eigenvalues 
A3 and A4. Note that the histogram is not symmetric with respect to the change Ai <H- A4 
and A® in -H- A® ax , since the eigenvalues are ordered, so the mean distance of the smallest 
eigenvalue to zero is smaller than the mean distance of the largest eigenvalue to unity. 
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C. Non Hermitian case and Multipartite operators 

The above analysis can be extended in a natural way for Hilbert spaces with m-fold tensor 
product structure, used to describe quantum systems consisting of m subsystems, 

U N = H ni ®---®Hn m , (30) 

with N = ri\ . . . n m . In the case of an operator X acting on this space, its product numerical 
range consists of all expectation values (Vw>d|^l?/Vod) among pure product states, \ip pm d) = 
\<f> m ). 

If the number m of subsystems is larger than two, there exist operators for which product 
numerical range forms a set which is not simply connected [TTJ, [22] . In fact the genus of this 
set can be greater than one. To show an illustrative example, we consider a unitary matrix 
of size two 

1 


The product numerical range of U® n can be found analytically for any integer n by applying 
an extension of the formula (|8]) to multipartite systems. Numerical range of U forms an 
interval / joining the complex eigenvalue e 1 ^ with the unity. Thus to find X®(U® n ), it 
suffices to compute the n-fold Minkowski power of the interval I on the complex plane. 
More explicitly, A®(?7® n ) consists of all the points z\z<i ■ ■ ■ z n , where z\ = 1 — A« + A^e"^ and 
Xi G [0, 1] for all i G {1,2,..., n}. Let us denote by / (a) the modulus \z\ of z = 1 — A + Ae l * 
as a function of the phase a := Arg(z). Obviously, / is a convex function of a. One can 
relatively easy get an explicit expression for /, 

cos (0/2) 

J K ' cos {a - 0/2) V ; 

Thus the numbers z of the form z = 1 — A + Ae*^, A G [0, 1] have a parametrization 



U 



(31) 



a i — y e ta f (a) with a G [0, 0] and / given by formula (32). Because of the convexity of /, for 
a fixed Arg (-21-22 . . . z n ), the minimum of \z1Z2 ■ ■ ■ z n \ is attained when z\ — z-i — . . . = z n . 
The resulting curve marks the border of A g, (f/ (X,n ) and has a parametrization 

\cos (a — 0/2) / 

The remaining parts of the border of A® {U® 11 ) are included in the n seg ments { [e^" 1 ^, e ik(i> ] }^ =r 
This follows because the maximum of \z\Z? ■ ■ ■ z n \ for a fixed (3 = Arg (z±z 2 . . . z n ) is attained 
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for Z\ = z 2 



z k -i = e**, z k = k<t>) and z k+l = z k+2 



z n — 1, where 



k — \_(3 /4>\. In Figure [4] we choose 4> = ^ and plot the product numerical ranges of U® n for 



n = 1,2,. ..,8. 




n = 5 



n = 6 



n = 7 



n = 8 



FIG. 4. Product numerical range of JJ® n for {7 specified in (31 ) with <j) = 3tt/5 and n = 1, . . . , 8 




FIG. 5. Product numerical range of U® 7 (plotted in gray) forms a set of genus 2. The matrix U 



is given by Eq. (31) with 4> = 3tt/5. 



Observe that if does not belong to the numerical range of U, it does not belong to 
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the product numerical range of U® n . Hence if for sufficiently large exponent n the product 
range 'wraps around' zero, the set h®{U® n ) is not simply connected. 

As one may notice in Figure |4| it is possible to construct a tensor product of operators such 
that its product numerical range has genus two. If we magnify picture number seven from 
Figure[4| it becomes evident that the genus of A®(f/® 7 ) is equal to two (cf. Figure[5]). Observe 
that if n is further increased, the genus of A®(U® n ) is not smaller than one, although the size 
of the hole around z = shrinks exponentially fast. More precisely, the distance between 
the set A (U) and zero is cos (0/2), which implies that the distance between A {U® n ) and 
zero equals [cos (0/2)] n for arbitrary n. 

In general, finding the product numerical range of a non-Hermitian operator without the 
tensor product structure is not a simple task. However, in the special case of a normal oper- 
ator X, which can be diagonalized by product of unitary matrices, a useful parameterization 
of its product numerical range was described in [22] . 



III. SEPARABLE NUMERICAL RANGE 



Consider a tensor product Hilbert space Hn = Hk ®Hm and the set f2 of all normalized 
states acting on it, p 6 Q p = p\ p > 0, Trp = 1. One distinguishes its subset f2 sep 
of separable states, i.e. states that can be represented as a convex combination of product 
states, 

pE Q scp ^ pen and p = '$2 Pi pf ] <8> p ( t I] (34) 

i 

Here positive coefficients pi form a probability vector, while pf^ and p[ denote arbitrary 
states acting on %k and 'Km-, respectively. Any state p which cannot be represented in the 
above form is called entangled [13] . Hence this definition depends on the particular choice 
of the tensor product structure, %n = Hk ® 'Km- 

Observe that Definition [6] of the product numerical range of an operator X acting on 
Wk ® "Hm can be formulated as 

A 9 (X) = {TrXp : p = \^ A ) <g> \^ B )^ A \ <g> ty B \} . (35) 

It is then natural to introduce an analogous definition of separable numerical range 

A scp (X):={TrXp:pen scp }. (36) 
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Since any product state is separable, the product numerical range forms a subset of the 
separable numerical range, A®(X) C A sep (X). By definition, the set fl sep of separable states 
is convex. This fact allows us to establish a simple relation between both sets. 

Proposition 1 Separable numerical range forms the convex hull of the product numerical 
range, 

A sep (X) = co(A®(X)) 

Proof. Assume that A G co(A®(X)), so it can be represented as a convex combination of 
points belonging to the product numerical range, A = YliPi^i- Taking the convex combina- 
tion of the corresponding product states = \ipf) <8> \4>f) we get a separable mixed state 
p = ^2iPi\(f>i)(4>i\ such that TrXp = A. A similar reasoning shows that if A ^ co(A®(A)) 
there is no separable state p such that TrXp = A. ■ 

Following [22] one can note that if A or B is normal then A sep (A (g) B) — A(A eg) B). 

Since product numerical range of a Hermitian operator forms an interval, in this case 
the separable and product numerical ranges do coincide. This is not the case in general. A 
typical example is shown in Figs. [6]d and ^jp, in which the separable numerical range forms 
a proper subset of the standard numerical range and includes the product numerical range 
as its proper subset. 

Consider, for instance a unitary matrix U of size 4 with a non-degenerate spectrum. 
Its numerical range is then formed by a quadrangle inscribed into the unit circle. If all 
eigenvectors of this matrix are entangled, the product numerical range of U does not contain 
any of its eigenvalues. In a generic case A® (£7) is not convex and it forms a proper subset of 
A sep ([7) - see Fig.§ 

A. k Entangled numerical range 

Any pure state in a iV = KM dimensional bipartite Hilbert space can be represented by 
its Schmidt decomposition, 

KM K 

m = EE^i*> ® \fi = Ev^i i ') ® io- (37) 
i=i j=i i=i 
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FIG. 6. Numerical range (light gray), separable numerical range (dark gray) and product nu- 
merical range (black dots obtained by random sampling) of family of matrices X a = C/^(a,0,0) • 



diag(i, — 1, — i, 1) • Ud{a, 0, 0)^, where U a is given by Eq. (26) for a = 0, tt/8, 37r/16, 7r/4. In the 
case of a = the eigenvectors of X form orthonormal canonical basis and X is normal therefore 
A scp (A) = A(X). In the case of a = it /A all eigenvectors of X are maximally entangled states and 
A se P(X) = A®(X). 



We the have assumed here that K < M and denoted a suitably rotated product basis by 
<g> The eigenvalues /ij of a positive matrix AA^ are called the Schmidt coefficients 
of the bipartite state The normalization condition \tp\ 2 = = 1 implies that 

II^IIhs = trAAf = 1, so the Schmidt coefficients \ii form a probability vector - see e.g. [IB] . 

The state |^>) is separable iff the K x M matrix of coefficients A is of rank one, so the 
corresponding vector of the Schmidt coefficients is pure. A given mixed state p is called 
separable if it can be represented as a convex combination of product pure states. This 
notion can be generalized in a natural way, and in the theory of quantum information [19] 
once considers set Vl^ of states which can be decomposed into a convex combination of 
states with the Schmidt number not larger than k. In symbols, p = Ylij=iVj\ ( t ) j){ ( t ) j\i with 
all vectors \<pj) = Yli=i ^il^Pf) ® H>f) °f Schmidt rank at most k. We may choose k to 
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be 1, . . . ,K, where K = M denotes the dimensionality of each subsystem. By definition, 
f^ 1 ) = Q scp represents the set of separable states, while = Q denotes the entire set of 
mixed quantum states. 

Making use of the definition of the subset of the set of all states in ^ one obtains 
an entire hierarchy of restricted numerical ranges denoted by A^. As the elements of 
are called /c-entangled states [2S], the set A^ k \X) will be referred to as numerical range 
restricted to fc-entangled states. 

For k = 1 one has fiW = ^ sep so in this case one obtains the separable numerical range, 
A^ 1 ) = A sep . Note that in this convention a 1-entangled state means a separable state. In 
the other limiting case k = K, Qk = & and one arrives at the standard numerical range, 
AW = A. The following chain of inclusions A® C A (1) C A (2) C • • • C A (K) = A holds 
by construction. This implies inequalities between the corresponding restricted numerical 
radii, r® < < <...,< = r. 



IV. APPLICATIONS IN QUANTUM INFORMATION THEORY 

In this section we link various problems in the theory of quantum information processing 
which have one thing in common: they can be analyzed using the restricted numerical range 
or related notions. 



A. Block positive matrices and entanglement witnesses 

Let us start by recalling the standard definition of block-positivity [13]. A Hermitian 
matrix X acting on the tensor product Hilbert space, Hn = Wm ® 7~t-K, is called block 
positive, if it is positive on all product states. Making use of the notation introduced in 
Sec. [iljthis property reads, A® in (X) > 0. Therefore checking if a given Hermitian matrix is 
block-positive is equivalent to showing that its product numerical range forms a subset of 
[0,oo). 

Block positive matrices arise in a characterization of positive quantum maps by the 
theorem of Jamiolkowski [3D] . A map $ taking operators on 1-L K to operators on %m, is 
called positive, if it maps positive operators to positive operators. Let denote 
the orthogonal projection onto the maximally entangled state \^?+) = -t= Yld=i acting 
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on Hk <8> 'H.k- The Jamiolkowski theorem states that $ is positive iff the corresponding 
dynamical matrix (Choi matrix [21]), Dq> = ($ ® 1) is block positive. This leads 

us to the following characterization of positive maps in terms of the product numerical range 
of _D<j>, 

Proposition 2 Let $ be a linear map taking operators on Hk to operators on Hm- Then 

$ is positive A®(£>$) C [0,oo). (38) 

That is, the product numerical range of Dq, has to be contained in the positive semiaxis in 
order for <3> to be positive. As discussed in Section [TTT] , for any Hermitian D, its product and 
separable numerical ranges do coincide. Consequently, positivity of $ can be formulated 
with A scp (/}$). The positivity condition reads: TrD^p > for any separable p. This is the 
same as A se P (£>$) C [0; +oo). 

We recall that a map $ is called fc-positive if $ £g> 1^ is a positive map. If this is the case 
for arbitrary k e N the map is called completely positive. The famous theorem by Choi [3T] 
concerning completely positive maps can be expressed in a similar manner, 

Proposition 3 Let $ be a linear map taking operators on Hk to operators on Hm- Then 

$ is completely positive <^ A (£>$) C [0, oo). (39) 



The difference is that (38) refers to the product numerical range of _D$ whereas (39) concerns 
the standard numerical range. Note that A (D$) C [0, oo) is just another way of writing 
that D$ is a positive operator. 

Positive maps find a direct application in the theory of quantum information due to 
a theorem by the Horodecki family [32]: a state a of a bipartite system is separable iff 
($ ® t)o~ > for any positive map $. In the opposite case, the state a is entangled. 

The above results explain recent interest in characterization of the set of positive maps. 
A block positive matrix W := Dq, which corresponds to a map which is positive but not 
completely positive, is called an entanglement witness, since it can be used to detect quantum 
entanglement. As discussed in sec. |III[ product and separable numerical ranges coincide for 
Hermitian operators. Thus the set of entanglement witnesses consists of Hermitian operators 
W such that TiW p > for all separable p and there exists an entangled state a such that 
Triy cr < 0. The set of separable quantum states can thus be characterized by a suitably 
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chosen set of entanglement witnesses. Such an approach was advocated in a recent work by 
Sperling and Vogel |33j, in which various methods for obtaining the minimal product value 
A® in of Hermitian matrices were analyzed. 

Bound [9] implies that the spectrum of an entanglement witness for any state of a K x K 
system has at most (K — l) 2 negative eigenvalues, in accordance with recent results of 
Sarbicki [31]. In the simplest case of K = 2, one recovers the known statement that any 
non-trivial entanglement witness in the two-qubit system has exactly one negative eigenvalue 



Our study of product numerical range of a Hermitian operator can thus be directly applied 
to the positivity problem. For instance, consider the family of one-qubit maps described by 



the dynamical matrix D = D(a,b,c) defined in (22). It is clear that these matrices are 
block positive iff G < 1/2. Therefore the expression (24) for G = G(a,b,c) gives us explicit 
constraints under which the map corresponding to D(a,b,c) is positive. If this map is not 
completely positive, the matrix D can be used as a witness of quantum entanglement. 

In the above case corresponding to maps acting on 2 dimensional Hilbert space H2 any 
2-positive map is completely positive. This is a consequence of the theorem of Choi [31], 
which implies that if a map acting on K dimensional Hilbert space is K positive, it is also 
completely positive. Thus for maps acting on a i^-dimensional system, it is interesting to 
study A;-positivity for k = 1 (equivalent to positivity), k — 2, . . . K — 1 and k = K (complete 
positivity). In general fc-block positive matrices are related to fc-positive maps. We are thus 
in a position to formulate the generalized Jamiolkowski-Choi theorem [2"5l making use 
of the concept of the restricted numerical range. 

Proposition 4 Let $ be a linear map taking operators on Hk to operators on Hm- Then 

$ is k-positive <^> A (fc) (£>$) C [0, 00). (40) 

As we explain in the next section, a special case of Proposition [4] for k = 2 is of relevance 
to the distillability problem for quantum states. 



B. n copy distilability of a quantum state 



It has been known for a long time (37] that bi-partite states with distillable entanglement 
are closely related to 2-positive maps and hence to 2-block positive operators (cf. also 
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[29, 38J). The precise relation between distillability and 2-block positivity is the following. 
Let p be an arbitrary state on a bipartite space Hn = Hk ® Wm- Assume that we allow 
only LOCC operations on a single copy of p. The state can be distilled into a maximally 
entangled state only if the partial transpose (1 <g> T) p is not a 2-block positive operator, 
i.e. it is not positive on states with Schmidt rank 2. Otherwise, p is one-copy undistillable. 



Writing this in terms of fc-entangled numerical ranges (cf. Section III A), we get the following 
proposition. 

Proposition 5 A state with a density matrix p on a bi-partite space Hn = Hk ® Hm is 
one-copy undistillable if and only if the 2-entangled numerical range of its partial transpose 
is contained in the nonnegative semiaxis, ((1 ® T) p) £ [0; +oo). 

If a state p turns out to be one-copy undistillable, it is still possible that a number of copies 
of p can be used for entanglement distillation. Proposition [5] is easily generalized to that 
situation. 

Proposition 6 Let p correspond to a state on a bi-partite space Hn = Hk ZSi'Hm- For any 
integer n, the state is n-copy undistillable if and only if the 2-entangled numerical range of 
(11 £g> T) p® n is contained in the nonnegative semiaxis, h.^ ((1 ® T) p® n ) £ [0; +oo). 



The symbol in Proposition prefers to positivity on states of Schmidt rank 2, where the 
Schmidt rank is calculated w.r.t. the splitting Hf/ 1 = 1-ifJ 1 £g> "Hf™ of the multipartite space. 
This is important to notice because many different splittings of T-L% n into a tensor product 
of two factors are possible. Evidently, Proposition [6] is nothing but Proposition [5] applied 
to p® n in place of p. This is easy to understand because the tensor product p® n represents 
a number n of identical, independent copies of the state p, e.g. coming from a source that 
produces p. 

It is natural to mention here a fundamental question concerning distillability of quantum 
states. Using the language of numerical ranges, we can formulate the problem in the following 

way: 

Given a density matrix p on a bi-partite Hilbert space H% n = U% n ®'H M n s.t. A ((1 <g> T) p) £" 
[0; +oo), can we infer that A^ 2 ^ ((11 £g> T) p® n ) ^ [0; +oo) for some positive integer nl 

In other words, is a bi-partite state p with a negative partial transpose always distillable, 
possibly using a huge number n of copies of p? This question has not yet been answered, 
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despite a considerable effort and some partial results (cf. e.g. [39|). 



C. Minimum output entropy and product numerical range 

Consider a completely positive map $ acting on the set fl^ oi normalized quantum states 
of dimension N. Minimum output entropy (see e.g. [4T)l Chapter 7]) is defined as 

S min ($) =min{5($(p))}, (41) 
p 

with p G Qn- Since the von Neumann entropy is concave, the minimum is attained on the 
boundary and thus 

5 nJn (*)= rmn{Sm^)(m}, (42) 

m(i>\ 

where G Qn are pure states. Therefore the minimum output entropy can be inter- 

preted as a certain measure of decoherence introduced by the channel. 

In (HI H2] it was proven that minimum output entropy (and thus Holevo capacity) is 
additive for unital channels. It is now known however that minimum output entropy is not 
additive in the general case (13] . Here we provide a characterization of the minimum output 
entropy for one-qubit channels using product numerical range of the dynamical matrix. 

Proposition 7 Let $ be a CP-TP map acting on f2 2 . Then 

S min ($) = A log(A) + (1 - A) log(l - A), (43) 

where X is a minimal value of product numerical range for the dynamical matrix Dq, 

A = A£ ta p,). (44) 

Proof. Let us define f(x) = —a; log 2 (x) — (1 — x) log 2 (l— x), which is increasing for x G [0, |]. 
Directly from the definition of minimum output entropy we can write 



S min ($) = mmS(<f>(\i)(i\)) = min/(A min ($(|z)(i|))) = / min A min ($(|i)(i|)) . (45) 

\V \i) V \V J 

Now since (k\Q(\i)(j\)\l) = (k <S> i\D$\l <S> j) (see [131 Eqn. 11.25]) we can rewrite the above 
expression as 

S^ n (®) = f(mm(j\<S>(\i)(i\)\j)) =/( f minO'®iP«|j®0 > ) = / (KJD*)) . (46) 
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Using the above proposition, we can easily calculate minimal output entropy for channels 
listed below. 

First we consider the amplitude damping, phase damping, phase flip, bit-flip and and 
bit-phase flip channels. In all those cases we can see from the Kraus form that the spectrum 
of the dynamical matrix has two zero eigenvalues. Then the plane spanned by the two 
eigenvectors corresponding to the zero eigenvalue contains at least one product state [25J. 
Thus A® in (.D<i>) = and the minimum output entropy for this channels is equal to zero. 
Using Proposition[7]one can also easily calculate minimum output entropy for some other one 
qubit channels. Consider the Werner-Holevo channel, described by the following dynamical 
matrix, 
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(47) 



for p £ [—1, 1/3]. In this case A^? n (D$ 

S min (*hw) = ~(1 - bl) log 2 |(1 - b|) - \{l + b|) log 2 1(1 + bl) (48) 

togfi-x) +2ptanh-» 

• (49) 



log(4) 

In the case of higher dimensional quantum channels we can use properties of product 
numerical range to check, whether for a given channel its minimal output entropy is equal 
zero. 

Proposition 8 For any completely positive, trace preserving (CP-TP) map we have 

S min m = iff 1GA%). (50) 
Proof. Since 1 G A® (/)$), there exists |z), \j) such that 

1 = (i®j\D*\i®j) = (51) 
Because $ is CP-TP channel, we have tr$(|j) = 1 and thus 

*(W01) = I*M (52) 
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The proposition follows. ■ 



D. Local discrimination of unitary operators 



The problem of local distinguishability of multipartite quantum states was analyzed by 
Walgate et al. [44J. Following their work, Duan et al. have shown [16] that two unitary 
operations U\ and U 2 are locally distinguishable iff G A®(V) where V = U\lJ 2 . If this is 
the case, then there exists a product state \ip) = \4>a,(Pb) such that the states Ui\ip) and 
U 2 \ip) are orthogonal and thus distinguishable. 

Our results on product numerical range allow us to solve the problem of local distin- 
guishability for a wide class of unitary operators. If the operator V = U\U 2 has the tensor 
product structure, V = V\ (g> V 2 , the two unitaries U\, U 2 are distinguishable iff the numerical 
range of any of the factors Vi, V 2 contains zero. This is the case when belongs to the convex 
hull of the spectrum of the factor V\ or of the factor V 2 - 




b) 




FIG. 7. a) Product numerical range for the matrix (53) with <f> = 2n/3 and ip = Wir/7, b) The 
region in the space of parameters ((p,ip) corresponding to locally distinguishable pairs (E/i,^) 



when U\U 2 =V(<f>,il>). 



Let us now deal with a more general case of V without the tensor product structure. 
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Consider for instance a family of unitary matrices of order four, 



10 

e** 



1 



(53) 



for (fi,ip E [0,2?r]. 

It is easy to show that the product numerical range of V is a bounded region of C whose 
border consists of the segments [e*^, l] , [l, e*^] and the line 



7 : [0, 1] 3 t ^ t 2 e llp + (1 - tf e ilp + 2t (1 - t) e C . 
For example, Fig. [T^i shows the shape of the product numerical range of V(^-, 



(54) 



107T ' 

7 - 



Using Eq. (54), it is not difficult to check for which values of the phases and if) 
the product numerical range of V (4>,if)) contains 0, so any U\ and U 2 such that U\U 2 = 



V ((f), if)) are locally distinguishable. Figure [7p) shows, in grey, the set of parameters (0, if)) 
corresponding to such distinguishable pairs (Ui,U 2 ). Explicitly, we have 



G A® (V (0,-0)) 1 1 sin if)\ cos0 + | sin0| cos-0 + 2y/\ sin^sin^l < 

A sin sin -0 ^ OA (0, if)) <£ {(0, 0) , (2vr, 2vr)} \. (55) 



For any two unitary matrices U\ and U 2 such that V = U\U 2 satisfies the above con- 
straints, it is possible to find a product state \x,Q with the property (x, £\V\x,£) = 0. A 
detailed construction of this state, presented in Appendix [Aj allows one to design the scheme 
of local discrimination between the unitary gates U\ and U 2 . 



E. Local fidelity and entanglement measures 

Several tasks of quantum information processing relay on the ability to approximate a 
given quantum state Qi by some other state g 2 . Alternatively, one attempts to distinguish 
Qi from q 2 . To characterize both problems quantitatively one may use fidelity, which can be 
interpreted as a 'transition probability' in the space of quantum states [45J, 



F{Qi,Q 2 ) = [Trlv^y^l)] 2 . 



(56) 
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We are going to follow here the original definition by Jozsa |46j, but one has to be warned 
that some later articles use the name 'fidelity' for y/~F. If one of the states is pure, Q\ = 



formula (56) simplifies and F = (V'll^lV'i)- Thus in this case fidelity has a simple 



interpretation of probability that the state Q2 is projected onto a pure state 

Consider two arbitrary mixed states Q\ and Q2 acting on a Hilbert space T-Ln- Although 



fidelity between these states is fixed and given by (56), one may pose a question to what 
extent fidelity can grow if local unitary operations are allowed. In other words, one asks 
about the fidelity between Q\ and U Q%U^ maximized over all unitaries U G U(N). This 
problem was studied in [47J , where the following bounds were established 

F(p\<r) < F( ei ,Ug 2 U j! ) < F(pV) = F(p\q^. (57) 

The vectors p and q represent the spectra of Q\ and Q2, while the up/down arrows indicate 
that the eigenvalues are put in the non-decreasing (non-increasing, resp.) order. Arguments 
of the fidelity in the above equation denote thus diagonal matrices which represent classical 
states. 

In this section we analyze an analogous problem for multipartite systems: What maxi- 
mum fidelity between two given states of such a system can be achieved, if arbitrary local 
unitary operations are allowed? We provide a solution of this problem in the special case 
when both quantum states are pure and derive bounds for the local fidelity in the case where 
p is a diagonal mixed state. 

Let |0) be a vector and g an arbitrary mixed state, both on Hab = Ha ® Hb- For 
simplicity we will restrict our attention to the symmetric case and assume that dim{^-L A ) = 
dim(n B ) = N. 

The fidelity of a mixed state with respect to a pure state is given by an expectation 
value, F = (4>\g\4>). We are going to study the question to what extend this quantity can be 
increased by applying arbitrary local unitary operations Ua®Ub- In other words, we look 
for the local fidelity defined as the maximum 

F max (f?,0) = max {<f>\{U A ®U B )^Q{U A ®U B )\(i>). (58) 

It is instructive to relate this quantity to a generalized numerical radius of an operator 
X, defined as the largest modulus of an element of its numerical range. Similarly for an 
operator X acting on a composed Hilbert space one defines product numerical radius as the 



27 



largest modulus of an element of A®(X). This notion can be further generalized, and for 
any operator X and an auxiliary operator C acting on the Hilbert space Hn — Ttk ® T~Lmi 
one defines the C-product numerical radius [T7] . 

r%(X) = max{\z\ : z = tr(E7i ® U 2 )X{U X ® £7 2 )t(7, ^ G E/(tf), t/ 2 e C/(M)}, (59) 

and other notions listed in Tab. |TTJ The problem of finding the local fidelity is then equivalent 
to determining the C-product numerical radius of the operator X = \<f>)((/>\ for C = g. 
Let us first solve the problem in the special case where the analyzed state is pure, g = 
It is then useful to represent both pure states using their Schmidt decompositions 



qsrp, 

N N 

10) = Yl v^> ® I*), i^) = Yl Vtoti) ® ( 6 °) 

t=l j=l 

The vector A of Schmidt coefficients set in a decreasing (increasing) order will be denoted 
by and A^, respectively. This notation allows one to formulate the following lemma. 

Lemma 1 For arbitrary local unitary operation Ua®Ub and pure states \<f>), one has 
< \(ip\U A ® U B \<fi)\ 2 < F(fj, 1 , \ l ) = I J Aj/ij ] . (61) 




The lower bound is a trivial consequence of the definition of fidelity. The upper bound 
follows from the theorem of Uhlmann which states that fidelity is given by the maximal 
overlap between purifications of both states, and the bound in [171 Eq. (4.19)]. This result 
follows also from the recent work of Schulte-Herbriiggen et al. [HI Prop. IV. 1]. 

If one of the states is separable, \<p) = \4>a) ® \4>b), its Schmidt vector has only a single 



non- vanishing component, A^ = (1,0, ... ,0), so the overlap (61) is bounded by the largest 
Schmidt coefficient /v ax of the state \ip). This is a special case of the geometric measure of 
entanglement of a multipartite state defined as the logarithm of the maximum projection 
on any product state [4"8] . 



E g (\7p)) = -\og(max\(7p\U l ®U 2 ®---®U m \Q,...a)\ 2 ) . (62) 

Here |0, . . . 0) represents an arbitrary product state, so transforming it by a local unitary 
matrix one explores the entire set of separable pure states of the m-partite system. Observe 
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that the argument of the logarithm in the above expression is just equal to the product 
numerical radius of the projector onto the analyzed state, X = 

In recent papers [IHJ ED] it was shown that the above maximization procedure becomes 
simpler if the multipartite state is symmetric with respect to permutations of the sub- 
systems, and all its coefficients in the product basis are non-negative. Then the maximum 
in (62) is achieved for the tensor product of a single unitary matrix, U\ oc = U® m , so the 
search for E g (%p) can be reduced to the space of a smaller dimension. It is then natural to 
ask whether this observation can be generalized for the problem of determining the product 
numerical radius of any multipartite Hermitian operator X, provided that X is symmetric 
with respect to permutations and it satisfies suitable positivity conditions. This problem 
was considered in a very recent paper by Hiibner et al. |51j . 

Thus the product numerical radius is useful in characterizing quantum entanglement of 
a pure state of a multipartite system. Interestingly, the product C-numerical radius of a 
Hermitian bi-concurrence matrix introduced by Badziag et al. [52] can be applied to describe 
the degree of quantum entanglement for any mixed state of a bipartite system. 

Let us then return to the bipartite problem and discuss the case when one of the two 



states in the expression (58) for local fidelity is pure while the other is mixed. Assume that 
the pure state \<p) is given by its Schmidt decomposition (60), while g is a diagonal mixed 



state, q = J2fj=i Pij\i) (*l ® The maximal local fidelity between these states can be 

bounded by the following lemma, proved in Appendix [Bj 

Lemma 2 The maximal fidelity between a pure state ip and diagonal state p is bounded from 
above, 

N 

- (63) 



max.)] pijB it 

ij=l 

where the maximum on the right-hand side is taken over all collections of non-negative real 



max F ((U® V)\i/>) (i/>\(U <g> V)\ g) < 

U,V&J{N) V ' 



numbers B^ that satisfy the constraints, for any {ii, «2 

N 

and for any {j 1 ,j 2 , ■ ■ ■ ,js} C {1,2,. 

N 



v}c{l,2,...,iV} 



■ + B irj ) 
N} 



N 



k=l 



k=N-r+l 



(64) 



Bij.] e 



N 



V(fc) 



fc=l 



k=N-s+l 



(65) 



where Am, A( 2 ), . . . , A(jv) are Schmidt coefficients of ip in ascending order. 
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The maximum on the right-hand side in Lemma [2] is attained at the edges of the polygon 



defined by the constraints (64) and (65). The bounds obtained in this way can be easily 
computed numerically using the simplex algorithm. 



F. Local dark spaces and error correction codes 

Consider a quantum operation $ acting in the space of mixed quantum states of size N, 
which can be represented in the Kraus form 

M 

p' = $(p) = ^r,pF/. (66) 
i=i 

To assure that the trace is preserved by the operation, the set of M Kraus operators has to 
satisfy an identity resolution, J^fLi Y}Yi = 1. 

Consider a /-dimensional subspace Pj = Yl\=i embedded in T-Ln- If it satisfies the 
set of M conditions 

P,X m p = A m P, for m = l,...,M (67) 

where X m = Y^Y m and X m G C no information goes outside of this subspace [55] , so p is 
called a dark subspace 



If a subspace P^ fulfils even stronger conditions of the type (67) 



PtfYjPi = X t ,Pi , for 2,j = l,...,M, (68) 

then quantum information stored in the system can be recovered, so the subspace Pi provides 
an error correction code USB] . Note that Pi has to simultaneously satisfy all the M 2 



equations (68). The complex numbers Ay corresponding to different X„-'s may be different. 



From an algebraic perspective condition (67) implies that A m belongs to the numerical 



range of order I of the operator X m [57]. In full analogy to the product numerical range, 
one may introduce the concept of product numerical range of higher rank as defined in Tab. 
II. This notion can be used to identify dark spaces or error correction codes with a local 



structure [58]. The distinguished subspace Pf, which solves the set of equations (66), can 
be chosen to be in the product form, P® = Yl\=i I* ® ® *!• 
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V. CONCLUDING REMARKS 

In this work we investigated basic properties of numerical range of an operator restricted 
to some class of quantum states. In particular, we analyzed the case of operators acting on 
a Hilbert space with a tensor product structure, often used to describe composed quantum 
systems. In this case one defines the product numerical range of an operator. We reviewed 
basic properties of this notion and presented some examples of operators for which product 
numerical range can be found analytically. 

To tackle the problem in a general case, however, we had to rely on numerical com- 
putations. In particular, we investigated an ensemble of TV = 4 random density matrices 
distributed according to the Hilbert-Schmidt measure and compared the probability dis- 
tributions of both edges of the product range with probability distributions for individual 
eigenvalues. 

In the case of a non-Hermitian operator its product numerical range forms a connected 
set in the complex plane. In general this set is not convex. The product numerical range of 
an operator acting on a two-fold tensor product is simply connected. However, this property 
does not hold for operators acting on a space with a larger number of subsystems. For any 
operator with a tensor product structure its product range is equal to the Minkowski product 
of numerical ranges of all factors. The theory of the Minkowski product of various sets in 
the complex plane, recently developed by Farouki et al. [23], can thus be directly applied to 
characterize the product numerical range of operators of the tensor product form. In this 
way we managed to establish product numerical range of a unitary product matrix JJ® n . 

Numerical range can also be generalized by taking other restrictions on the set of quantum 
states. Although we studied here the case of numerical range restricted to separable and 
fc-entangled states, one may also use other restricted sets of quantum states or combine 
these conditions, analyzing for instance the set of real product states. As the product states 
of the K x M system can also be considered as coherent states with respect to the composite 
group SU(K) £g> SU(M) [59], an analogous relation holds for the corresponding numerical 
ranges. 

Numerical range can also be generalized in other direction: for each case of a restricted 
numerical range one can introduce concepts and generalizations known for the standard 
numerical range. In Table [IT] we have collected standard definitions of numerical range, 
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Standard definitions 
(for simple systems) 
X : Hn T~Ln 


Product definitions 
(for multipartite systems) 
X : n ni ® ■ ■ ■ ® H nm -»• H ni ® • • • ® ^n m 


numerical range 
A(X) = {( ( p\X\ ( p):\ ( p)eH} 


product numerical range 
A®(X) = {(V>i 8) . . . ® ^ m |X|^i ® ... (8) -0 m ) : |^) G WnJ 


numerical radius 
r(X) = maxjlzl : z G A(X)\ 

V / L 1 1 V / J 


product numerical radius 
r 9 (X) = max{|z| : 2 G A®(X)j 

V / L 1 1 \ / J 


C-numerical range 
A C (X) = {A : A = tiUXWC} 
where U G U(N). 


product C-numerical range 
A®(X) = |tr(C/i ® . . . ® U k )X(U! ® . . . ® CC+C} 
where C/j G U(rii). 


C-numerical radius 
r c {X) = max{|z| : z G A C (X)} 


product C-numerical radius 
rg(X) = max{|z| : z G Ag(A)} 


higher rank numerical range 
A t (X) = {A : PiXPi = XP t } 


higher rank product numerical range 

Af(X) = {A : PfXFf = AP®} 
i^ 8 = ELi 1 1® ; " ® i) (»<» •■•<»» 1 



TABLE II. Standard algebraic definition of the numerical range and related concepts compared 
with their product analogues. The definitions on the left concern an operator X acting on Hilbert 
space TIn, while their product analogues are defined for operators acting on a tensor product 
Hilbert space % ni ® • • • ® 7i nrn ■ Here \tp) denotes an arbitrary state of "Hat, while \ip\ ® • • • ® tp m ) = 
IV'i) ® ' ' ' ® |V>m) G 7ii ® . . . % m represents an arbitrary product state of the composite, m-particle 
system. 

numerical radius, C-numerical range and higher rank numerical range [1)7], along with their 
counterparts denned for Hilbert space of the form of an m-fold tensor product, Hn = 
H ni <2) • • • <8) Hnmi with N = ni . . . n m . Note that C-numerical range, as well as product 
C-numerical range, reduce to the numerical range (product numerical range, resp.) for 
C = diag({l, 0, . . . , 0}) and this case was already analyzed in [T5] . 

Observe that the above concepts arise naturally in a variety of problems in quantum 
information theory. For instance, being in a position to find the product numerical range of 
an arbitrary operator, one could advance fundamental problems concerning the characteri- 
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zation of the set of positive maps or description of the set of entangled states and finding 
the minimum output entropy of a one-qubit quantum channel. Therefore, improving tech- 
niques of finding restricted numerical ranges would have direct implications for the theory 
of quantum information and quantum control. For example, in this work we have estab- 
lished the positivity of a certain family of one-qubit maps, we solved the problem of local 
distinguishability between a class of two-qubit unitary gates and analyzed the properties of 
local fidelity between quantum states. 

In conclusion, we advocate further studies on restricted numerical range and cognate 
concepts. On one hand, the restricted numerical range is an interesting subject for mathe- 
matical investigations. On the other hand, it proves to be a versatile algebraic tool, useful 
in tackling various problems of quantum theory. 
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Appendix A: Product vectors for local discrimination between U\ and U2 



The discussion in Section [IV D| left aside the question of precisely how the unitaries U\, U2 
fulfilling A® (u{U 2 ^ = can be distinguished. To accomplish this task in practice, one needs 
to find a product vector \x,Q such that (x, Cl^i^lx, = 0- There exists in general an 
infinite number of such vectors. In the case U\U?, — V (0, ip) analyzed in Section 



IV D 



it is 



not difficult to find all of them. Recall that V (0, ■0) is of the diagonal form diag (l, e 1 ^, e 1 ^, l) 
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with respect to the tensor product basis (|00), 1 01) , |10), |11)} of H2 <8> 'H.i- Let us write 
|f) = y/ie iKo \0) + s/l^te iKl \l) and \x) = ^/se iSo \0) + s/l^se i5l \l) for s,t G [0,1] and 
k , Ki,5o,5i arbitrary real numbers. Thus we assume that \x) and |£) are of unit norm, 
which is permissible. It is now easy to see that 

(X, Z\V (0, V) \X, () = ts + (l- t) (l-s) + e«t (1-8) + (1 - t) s, (Al) 
where s, t G [0, 1]. 



Formula (Al) gives us some idea of how the results presented in Section IV D were ob- 
tained. Note that the phases k , ki,S and Si are irrelevant to the value of (x, £\V (0, ip) \x, £)• 
Thus any product vector that fulfils certain relations between the amplitudes t and s can 
be used for perfect discrimination between the two unitaries. Note that this is a general 
property whenever \j\Ui is diagonal with respect to some tensor product basis of %k ®Hm 
and G ^{u\uX 

In order to solve Eq. (Al) for s and t, we first observe that 

Im(< X ,£|V(0,V>)lx,£)) = O 

reduces to sin0t (1 — s) + sin^s (1 — t) — or 

t sin0 



(A2) 



t (sin + sin ip) — sin if) 



(A3) 



If we substitute this in Eq. ( Al ), we get the condition (x, £\V ((f), ip) |%, £) = in the following 
form 

(1 -t) (tsia((f)-ip) - (1 -t)siuip) = 0. 



t 2 sin ( 



(A4) 



We can solve (A4) for t G [0,1] under the assumption that G A® (V (0, -0)) (cf. the 



conditions on the right-hand side of (55)). The result is 



sin (0 — ip) 2 + 4 sin sin ip + | sin (0 — ip) \ + 2| sin -0| 
2 (| sin0| + I sin (0 — -0) I + I sin-0|) 
By symmetry we obtain an expression for s, 



(A5) 



sin (ip — 0) 2 + 4 sin if) sin + I sin (-0 — 0) I + 21 sin 0| 

(A6) 

2 (\sinif)\ + \sm (if) - (f))\ + \sin(f)\) ' K ' 

Hence the product vector useful for perfect local discrimination between Ui and U2 can be 

any of the family 



(yie iKO \0) + VT^te iK1 \l}^j g> (yGe^lO) + VT^e l5l \l}) 



(A7) 
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with Ko, Ki, So, Si G R and s,t given by the formulas (A6) and (A5), respectively. This only 
works when U\U 2 = V 0, ip) and G A® (u\U 2 



Appendix B: Proof of Lemma [2] 

Let us introduce matrix A which depends on the vector A and a local unitary matrix 
U <8> V, with entries 



Aij 



N 



j2Vh(u(k) \i)(v(k) \j) 



k=i 



(Bl) 



where by U (k) we mean U (\k)), so that (U (k) \i) corresponds to {k\U^\i) in the usual 
physicists' notation. Similarly, (k\U (i)) = (k\U\i). 



Using the notation of eq. (Bl ), we arrive at a handy expression for the expectation value 



N 



mu®vy P (u®v)\i/>) = Y / p ij A ij , 

ij=l 



(B2) 



which we wish to maximize over the set of local unitaries. The first thing to notice is that 
Aij are non-negative real numbers and 



A' 



(B3) 



— i, 

ij=l 

thus the matrix A treated as vector is an element of standard (N 2 — l)-simplex. 
Matrix A defined above satisfies the following lemma. 



Lemma 3 For any {i\,i2, ■ ■ ■ , v} C {1,2,..., N} 

N 

+ A i2 j + • • • + A irj ] G [A(i) + A( 2 ) H h A( r ), A(jv) + A(jv-i) H h A(Ar_ r+ i)] (B4) 

and /or any j 2 , . . . ,j s } C {1,2,..., A} 

N 

/ J [Aij! + A ij2 + • • • + A ijs ] g [A(i) + A( 2 ) H h A( s ), A(jv) + A(jv-i) H h A(jv- s +i)], (B5) 



where A(i), A( 2 ), . . . , Any) o^e tae Schmidt coefficients of \^) sorted ascendingly. 
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Proof. First we write 

N 

/y [Ahj + A i2 j + • • • + A irj ] = 

N [ N N 

= E i E ^(umiv (k) b)i 2 +---+iE y^( u (*) ( k ) \M 

j=l L k=l k=l 
N N 

= E[ E V^(U(ki)\ii)(V(k 1 )\j)(i 1 \U(l 1 ))(j\V(l 1 )) + 

j=l kih=l 
N 

■■■+ E V^K(U (K) \i r )(V (K) \j)(i r \U (l r ))(j\V (l r )) 
N N 

= E V^(u(ki)\ii)(ii\u(i 1 ))J2(v(h)\j)(j\v(ii)) + 

kxh=l j=l 

N N 

••■+ E V / VA^(f/(^)|^)(v|^(/r))^(^(A; r )|j)(j|y(/ r )). 

fc r / r =l j=l 

Since |j) form a basis, we have the identity ^2 \ — 1- Thus 

N 

E + ^ i2j ' + " ■ + ^vj] = 

AT 

= E V / ^A(^(^i)Ki)(n|f/(/i))(V(A;i)|V r (/:)) + 



fci/i=i 



iV 



■■+ E V^k r K(U (K) \ir)(ir\U (l r ))(V (k r ) \V (l r )) 



N 



N 



(B6) 



= K\(u(h) Ni)| 2 + ■ ■ ■ + E X ^\{U (Ar) |v)| 2 - (B7) 

fcl=l k r = l 

This fact, combined with Corollary 4.3.18 of Horn and Johnson [60J, proves the lemma. ■ 
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